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We systematically calculate the total reaction cross sections of oxygen isotopes, ^^"^''O, on a 
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determinants generated from a phenomenological mean-field potential. The agreement between 
theory and experiment is generally good, but a sharp increase of the reaction cross sections from '^^O 
to remains unresolved. To examine the sensitivity of the diffraction pattern of elastic scattering 
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at the incident energy of 300 MeV by calculating the full Glauber amplitude. 
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I. INTRODUCTION 

Studies on neutron-rich unstable nuclei have been at- 
tracting much attention both experimentally and theo- 
retically. These studies are motivated, for example, by 
that we want to understand the nuclear structure and ex- 
citation mode of the neutron-rich nuclei as well as a role 
played by them in forming heavy elements in stars. Bind- 
ing energy, radius and density distribution, among oth- 
ers, are basic quantities to determine the nuclear prop- 
erty. Reactions of unstable neutron-rich nuclei with a 
proton target are, therefore, of current interest as they 
are at present a major means to probe the matter densi- 
ties of exotic nuclei, particularly the region of the nuclear 
surface. If one appropriately selects the incident energies, 
protons could be more sensitive to the neutron distribu- 
tions than to the proton distributions of nuclei. 

The Glauber theory Q offers a powerful and handy 
framework for the description of high energy nuclear re- 
actions. This theory describes proton-nucleus reactions 
very well. Usually a calculation based on the Glauber 
theory is performed using a one-body density calculated 
from a nuclear wave function. However, Bassel et al. 
calculated Glauber's scattering amplitude using a Slater 
determinant wave function, and were able to examine the 
effect of Pauli blocking, which is impossible to discuss us- 
ing just the one-body density. 

Recently, we have analyzed the total reaction cross sec- 
tions of carbon isotopes on both ^^C Q and proton Q 
targets using the Glauber theory. The densities of the 
carbon isotopes are constructed from Slater determinants 
generated from a phenomenological mean-field potential 
or two types of dynamical models, core-|-n and core-|-n-|-n 
models, which go beyond the mean-field model. 

Data on the cross sections of the oxygen isotopes up 
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to the dripline nucleus '^^O are available at high ener- 
gies Q : The cross sections show a gentle increase with 
increasing neutron number from ^^O to ^^0 and then an 
abrupt increase up to ^"^O. On the other hand, several 
theoretical studies on the matter radii predict a mild in- 
crease of the radii with increasing neutron number. The 
calculations have been performed for the isotopes up to 
iV = 12 0, and much further 0, i, 0, [HI- As the 
spin-parity systematics of the ground states of the oxy- 
gen isotopes indicates that the single-particle orbit 
seems fairly stable, we first analyze the cross section data 
using a model similar to that of the previous study, that 
is, assuming the mean-field potentials that reproduce the 
nucleon separation energies. The interest of this analysis 
is to examine the extent to which the theory can explain 
the characteristic behavior of the cross sections. 

One of the advantages of the Glauber theory is that the 
same input as used in the reaction cross section calcula- 
tion is readily applicable for the calculation of the differ- 
ential cross section of clastic-scattering. The diffraction 
pattern in the differential cross section is expected to de- 
pend on the diffuseness of the nuclear surface. Therefore 
we also study elastic scattering of the oxygen isotopes on 
a proton target. The case of p-*^0 scattering serves as a 
testing ground of our model. We then predict the differ- 
ential cross sections of p-20'2i,23Q gja,stic scattering. The 
differential elastic-scattering cross section data of p-^°0 
is available only at low energy, such as ^ 30 MeV per 
nucleon p^ . 

This paper is organized as follows. The reaction mod- 
els for the calculation of the total reaction and elastic- 
scattering cross sections are recapitulated in Sect. [Ill The 
construction of the densities of the oxygen isotopes is ex- 
plained in Sect. IIIII We present the calculation of the 
total reaction cross sections of the oxygen isotopes in 
Sect. IIVI The differential cross sections for neutron-rich 
oxygen isotopes are given in Sect.|Vl A summary is given 
in Sect. I VII Appendix A presents the densities of the oxy- 
gen isotopes which are obtained using the Slater deter- 
minants constructed from the harmonic-oscillator single- 
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particle orbits. Appendix B explains the calculation of 
the p -nucleus elastic-scattering amplitude using a Slater 
determinant. 



II. PHASE SHIFT FUNCTIONS IN THE 
GLAUBER MODEL 



The Glauber theory provides us with an excellent 
framework to describe high energy reactions [l|. The 
p -nucleus elastic-scattering amplitude is given by 



(1) 



where b is the impact parameter vector perpendicular 
to the beam (z) direction, K the initial momentum of 
the relative motion, and q the momentum transferred 
from the projectile to the target. Here we neglect the 
p -nucleus Coulomb potential, because its effect appears 
only in the extreme forward direction for a light target. 
The nuclear part of the optical phase-shift function is 
given by 
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where ^0 is the intrinsic (translation-invariant) A- 
nucleon wave function of the nuclear ground state, and Si 
is the projection onto the xy-plane of the nucleon coor- 
dinate relative to the center-of-mass of the nucleus. Here 
is 1 for neutron and —1 for proton. 
The profile function, V^^^ for the and pn elastic 
scatterings, is usually parameterized in the form 



(3) 



where apN is the ratio of the real to the imaginary part 
of the pN scattering amplitude in the forward direction, 
a^°^ is the pN total cross sections, and (ipN is the slope 



parameter of the pN differential cross section. These 
parameters for different incident energies are tabulated 
in Refs. I,!]. 

The differential cross section of the p-nucleus elastic- 
scattering is given by 



da 
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\F{q)\ 
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and the total reaction cross section of the p -nucleus col- 
lision is calculated from 



db 1 



(5) 



In the optical limit approximation (OLA), the p- 
nucleus scattering phase-shift function is simply given 



using the proton density Pp{r) and the neutron density 
/9„(r) as follows 



aiXOLA(b) 



cxp [iXpib) + iXn(fe)] 
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where Xp iXn) denotes the phase shift due to the protons 
(neutrons) inside the nucleus 

iXpib) = - J drpp{r)Tpp{b- s), 

iXn{b) = - j drpn{r)Tpn{b - s). (7) 

We will evaluate Eq. ^ completely for a 'I'o that 
is given by a Slater determinant comprising harmonic- 
oscillator single-particle orbits, and use it in Sect. IVl 

III. DENSITIES OF OXYGEN ISOTOPES 

Compared to the carbon isotopes, where the compe- 
tition of the lsi/2 and 0^5/2 orbits appears to play an 
important role to determine their structure, the spin- 
parity systematics of the oxygen isotopes indicates that 
the Qd^/2 orbit is lower than the lsi/2 orbit in the ground 
states up to ^^O. The configurations for the oxygen iso- 
topes are thus assumed to be given according to the shell- 
model as follows: The nucleus -'^^O is a doubly magic nu- 
cleus occupying the Osi/2, OP3/2 and Opi/2 orbits, and 
^^O has a neutron hole in the Opi/2 orbit. The nucleus 
with 16< A<22 has A—\% valence neutrons in the 
0^5/2 orbit, and for A=23, 24 it has six neutrons in the 
0^5/2 orbit as well as ^ — 22 neutrons in the lsj^/2 orbit. 
We assume the spin-parity of the ground state of ^^O is 
l"*^ though it is not yet confirmed experimentally 

The single-particle orbits are generated from the mean- 
field potential containing central and spin-orbit poten- 
tials. The Coulomb potential is added for protons. The 
strength of the spin-orbit potential is set to follow the 
standard value, whereas the strength of the central part 
is chosen so as to reproduce the separation energy of the 
last nucleon. See Ref. Q for details. This prescription is 
called an "S'„ model" hereafter. This model apparently 
ignores the pairing effect, which gives larger separation 
energy for the even- TV nucleus than for the odd- TV nu- 
cleus. Therefore the S'„ model tends to predict a too 
large size for the odd-A^ isotope. To remedy this prob- 
lem, we also test another one, called an "(S'„) model", 
which fits the average separation energy for the nucle- 
ons in the last orbit. For example, in the case of ^^O in 
which three neutrons occupy the 0^5/2 orbit, the average 
neutron separation energy to be fitted is one third of a 
sum of the neutron separation energies of ^'^''iS'iSO, and 
likewise the average proton separation energy to be fitted 
is one half of a sum of the proton separation energies of 
^^O and ^^N. The center-of-mass motion is taken into ac- 
count in order to obtain the intrinsic density. See Ref. Q 
for detail. 
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FIG. 1: The matter, neutron and proton root-mean-square radii of the oxygen isotopes calculated with the Sn and (Sn) models. 
The empirical values are taken from Ref. 



We also use harmonic-oscillator single-particle orbits 
in Sect. |V]to examine the differential cross section. 

Figure [T] compares the matter, neutron and proton 
root-mean-square (rms) radii of the oxygen isotopes (r,„ , 
r„, and r^, respectively) calculated in the Sn and (5,1) 
models. Clearly the strong even-odd staggering seen in 
the neutron radius of the 5„ model becomes mild in the 
{Sn) model, and the matter radii in the latter model seem 
to be in better agreement with those extracted from a 
model-dependent analysis of the interaction cross section 
data Noticeable discrepancies appear at ^^O and ^^O. 
The theory predicts a smaller radius for ^^O but a larger 
radius for ^'^O than the empirical values. We will see 
in Sect. IIVI that this discrepancy directly appears in the 
comparison of the reaction cross sections. The charge 
radii of le.iT.isg calculated in the (5„) model using the 
finite size correction of the nuclcon are found to be 2.64, 
2.63, 2.62 fm, which are slightly smaller than the mea- 
sured values, 2.718(21), 2.662(26), 2.727(20) fm [l|, re- 
spectively. 



IV. RESULTS FOR TOTAL REACTION CROSS 
SECTIONS 

Before showing results of calculation for the total reac- 
tion cross sections of the oxygen isotopes, we first com- 



pare in Fig. [2] the theoretical total reaction cross 

section to experiment. The OLA approximation is used. 
The cross sections calculated at higher energies are in 
reasonable agreement with the data. A slight underes- 
timation is consistent with the fact that the density of 
1^0 gives a slightly smaller radius compared to the ex- 
perimental value. The theory predicts about 15% larger 
cross sections at lower energies, which is very similar to 
the one found in the analysis of the p-^^G total reaction 
cross section The reason for this discrepancy is due 
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FIG. 2: Comparison of the p-f ^''O total reaction cross sections 
calculated in the OLA with experiment. The data are taken 
from Ref. 
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FIG. 3: (Color online) Total reaction cross sections for the 
oxygen isotopes incident on a ^^C target. The incident energy 
is about 1 GeV per nucleon except for ^^O. The calculation 
is based on the NTG approximation. The experimental data 
are taken from Ref. [^. 



to the limitation of applying the Glauber theory to the 
reaction at the energy lower than 100 Me V [l6l |. 

Figure[3]displays the total reaction cross sections of the 
oxygen isotopes incident on a ^^C target. The calculated 
values are based on the NTG approximation which 
includes higher-order corrections missing in the OLA cal- 
culation, and they are compared to the interaction cross 
sections measured at incident energies around 1 GeV per 
nucleon 0]. The energy per nucleon is, however, about 
710 MeV in the case of ^^O. The density used for is 
the same as in Ref. The (S'„) model seems to give 
slightly better results than the Sn model. 

For a comparison between theory and experiment, it 
is important to realize that it is the total reaction cross 
section that is calculated and not the interaction cross 
section. The interaction cross section does not include 
the contribution from the inelastic excitation of a projec- 
tile nucleus to a particle-bound state, whereas the total 
reaction cross section includes this contribution. It is in 
general not trivial to estimate the difference between the 
two cross sections [l3|, but they are not expected to be 
very different at high energies [H, [ll| . As already noted 
in the comparison of the matter radii, the theory can 
reproduce the increase of the cross sections up to 
but fails to reproduce the sharp rise of the cross section 
from to 2^0. It was noted in Ref. ^ that the latter 
discrepancy cannot be explained in any theoretical mod- 
els. See Fig. 4 of Ref. [201 ^ '^sll as the calculation of 
Ref. 

The sharp rise in the interaction cross section at ^'^O 
has led to some controversial issues for the understanding 
of its origin. It was conjectured in Refs. [lO, [131 that a 
core nucleus ^^O could be substantially modified in ^■^O 
if described with a model and also that the spin- 

parity of the ground state of ^^O might be I"*" rather than 



. Motivated by this, several experiments have been 

performed, which aU support [H, [13, [H, HI, [13 , as 
adopted in the present study. Namely, the configuration 
of the ground state of ^^O is dominated by a lsi/2 neu- 
tron coupled to ^^O{0~^) and an excited configuration of 
(151/2)^(0^5/2)^^ appears as a resonance lying very close 
to the n+^^O threshold [H, [13|. The enhancement of 
the cross section from ^^O to ^'^O thus remains an open 
question to be explained. 

As another possibility for getting information on nu- 
clear size, particularly near its surface, we will examine 
in the next section the p -nucleus elastic scattering. 



V. ELASTIC-SCATTERING CROSS SECTIONS 

In the previous section we calculated the total reac- 
tion cross section for the oxygen isotopes using the phe- 
nomenological mean-field potential. Though the (Sn) 
model was found to give nuclear sizes that reproduce 
the cross sections reasonably well, the increase of the 
cross section from ^^O to ^^O was not explained. The 
p-2°0 elastic scattering has recently been measured at 
the National Institute of Radiological Sciences (NIRS) in 
Japan [1^ . In this section we study the differential cross 
section of the p-oxygen isotope elastic-scattering in order 
to examine how sensitive the differential cross section is 
to the nuclear size. 

To this end, it is more convenient if we can adjust 
the nuclear size more flexibly. Thus we use harmonic- 
oscillator single-particle wave functions instead of the 
orbits determined from the separation energies as was 
done in Sect. IIIII The size parameters Vp and f„ of the 
harmonic-oscillator wave functions are first determined 
so as to reproduce the proton and neutron rms radii of 
the (Sn) model and then, by changing the f„ value to 
fit the total reaction cross section, we examine the ex- 
tent to which the differential cross section is altered. Ta- 
ble [J lists the radii and the total reaction cross sections 
of the oxygen isotopes for several sets of i^p and Vn val- 
ues adopted in this section. The r^, r„, and rp values 
in each row containing a parenthesis correspond to the 
values of the (Sn) model. The density formula obtained 
from the harmonic-oscillator functions are summarized 
in Appendix A. 

The calculation of the Glauber amplitude ([1]) for the 
p -nucleus elastic scattering is explained in Appendix B. 



A. Test of the Glauber amplitude: "O 

We present numerical results of the differential cross 
section for the elastic-scattering, and compare it 

to experiment to study how precisely we can discuss the 
nuclear size within the Glauber theory. 

We start with the calculation of the differential cross 
section for p- 1^0 at an incident energy of 1.0 GeV, where 
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FIG. 4: (Color online) Comparison of the numerical results 
with the differential elastic-scattering cross section data for p - 
^^O at -Bp = 1.0 GeV. The solid curve represents the full calcu- 
lation with the harmonic-oscillator shell-model wave function. 
The dashed curve represents the calculation with OLA. The 
data are taken from Ref. [^. 
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FIG. 5: (Color online) Comparison of the numerical results 
with the differential elastic-scattering cross section data for 
at £p=317MeV. See the caption of Fig. H 



the parameters of the nuclcon-nucleon scattering ampli- 
tude are well determined. The numerical results of the 
differential cross sections are shown in Fig. [D The solid 
curve represents the result of the full calculation with 
the harmonic-oscillator shell-model wave function, while 
the dashed curve shows the cross section calculated in 
the OLA. Both results agree satisfactorily with the data. 
The two calculations give no large difference up to the 
second minimum. This analysis confirms that the 
elastic-scattering cross section is reproduced using the 
wave function which has a correct size together with the 
profile function used in the present study. 

As the experimental study of p -nucleus elastic scatter- 
ing at around 300 MeV is considered an interesting and 
promising project at RIKEN, it is useful to assess the 
suitability of the profile function at this energy as well. 
For this purpose we calculate the differential cross section 
of the elastic scattering in this energy region. Fig- 

ure O displays the results of the differential cross section. 
Both the exact and OLA curves reproduce the data fairly 
well, similarly to the case of 1 GeV. From these compar- 



isons at 300 MeV and 1 GeV, we can conclude that the 
p -nucleus elastic scattering can reliably be described in 
the OLA at least up to the second minimum using the 
present parameters of the profile function. 



B. The cross sections of ^i-^^o 

As was shown in the previous section, the total reac- 
tion cross sections of the oxygen isotopes with a ^^C tar- 
get are reproduced reasonably well, except for the cases 
of ^^'^•^O. In particular, the total reaction cross section 
of ^^0-1- "'^^C is calculated to be about 2% smaller than 
the measured interaction cross section. This may sug- 
gest that the (Sn) model slightly underestimates the mat- 
ter radius of ^^O. As mentioned at the beginning of the 
present section, we increase the radius to reproduce the 
interaction cross section. Three sets of Vp and h'n values 
arc listed in TablcIH In the case of ^"'^O, the total reaction 
cross section calculated with the NTG approximation is 
about 5% larger than the interaction cross section. As 
^^O has four particle-bound excited states, the difference 
between the two cross sections could be reduced to some 
extent. Table |T] lists a set which fits the interaction cross 
section of ^lO+^^C. 

Now we examine whether or not the p -nucleus elastic 
scattering can give us useful information on the matter 
radius of the nucleus. Figure [5] displays the differential 
cross section of elastic-scattering at i?p=300MeV 

for the three different matter radii listed in Table [D We 
focus on the angle of the first peak and the magnitude 
of the cross section at this angle. As the matter radius 
increases, the first peak moves to a smaller angle. The 
angle in degrees (the magnitude of the cross section in 
mb/sr) of the first peak for the solid, dashed, and dotted 
curves are 21.5 (4.7), 21.1 (4.2), and 20.5 degree (3.3 



TABLE I: The rms radii, in fm, of the matter, neutron and 
proton density distributions for several oxygen isotopes. The 
Vp and Vn values given in fm are the size parameters of the 
harmonic-oscillator functions for proton and neutron. See Ap- 
pendix A. The seventh and eighth columns denote the total 
reaction cross sections, in units of mb, calculated in the OLA 
for proton and '^^C targets at the incident energy of 1.0 GeV 
per nucleon. The interaction cross sections taken from Ref. 
are given in parentheses. 
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FIG. 6: (Color online) Comparison of the differential elastic- 
scattering cross sections of p-^^O at i5p=300MeV for differ- 
ent values of nuclear matter radii. The solid curve represents 
the case with a radius of 3.05 fm, the dashed curve with a 
radius of 3.15 fm, and the dotted curve with a radius of 3.25 
fm. The calculations are carried out with the OLA. 



mb/sr), respectively. 

The difference, at each first peak angle, between the 
magnitudes of the cross sections denoted by the solid and 
dotted curves is about 30%. This value is not small, but, 
unfortunately the presently estimated experimental un- 
certainty at this angle for the p-^°0 elastic scattering is 
comparable to this value [l^ . The angle shift of the first 
peak position between the two curves is about 1.0 degree, 
which might be too small to distinguish experimentally. 

We have carried out a similar analysis for the case of 
^^O. The result is displayed in Fig. [T] The angle (the 
magnitude of the cross section) of the first peak for the 
solid and dashed curves arc 23.5 (7.3), 22.3 degree (6.4 
mb/sr), respectively. Again a high precision experiment 
will be needed to distinguish which of the above two is 
favorable. 



C. Prediction of the cross section of '^"O 

Recently the differential cross section of ^"O elastic 
scattering on a proton has been measured at the energy 
of 300 MeV per nucleon at NIRS Motivated by this 
experiment, we predict the elastic and total reaction cross 
sections of p-^"0. 

The Glauber amplitude is calculated fully using the 
Slater determinant consisting of the harmonic-oscillator 
single-particle orbits. The size parameter of the harmonic 
oscillator wave function is chosen to reproduce the radius 
of ^"O given in the (Sn) model, as listed in Table HI The 
result is shown in Fig. [5] The first peak appears at 24.8° 
and the cross section at the peak is 2.43 mb/sr. As shown 



in Fig.[3l the reaction cross section calculated in the NTG 
approximation is about 30 mb larger than the interaction 
cross section. It is not clear, however, whether this dif- 
ference between the two cross sections indicates that the 
radius used for is slightly too large or not, because 
a considerable contribution to the reaction cross section 
is expected from the 12 particle-bound excited states of 
^"O. In any case the predicted differential cross section 
near the peak appears consistent with the preliminary 
data m. 



VI. SUMMARY 

We have systematically calculated the total reaction 
cross sections of the oxygen isotopes, "'^^"^^O, on a ^^C 
target at high energies using the Glauber theory. The 
oxygen isotopes are described with the Slater determi- 
nants generated from a phenomenological mean-field po- 
tential, which is an extension of our previous work for 
describing the carbon isotopes [3, We have intro- 
duced two schemes, the 5,1 model and the (Sn) model. 
The agreement between theory and experiment is gen- 
erally good, especially in the model, but the sharp 
increase of the reaction cross sections from ^^O to ^'^O 
remains an open question. 

As a possible cross section which may depend on the 
nuclear radius more sensitively than the reaction cross 
section, we have examined the differential cross sections 
of the p-20'^i'23Q elastic-scatterings at the incident en- 
ergy of Ep = 300 MeV using the full and OLA Glauber 
amplitudes. The differential cross sections calculated 
from the two amplitudes are not very different up to the 
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FIG. 7: (Color online) Comparison of the differential elastic- 
scattering cross sections of at i5p=300MeV for differ- 
ent values of nuclear matter radii. The solid curve represents 
the case with a radius of 2.64 fm, while the dashed curve with 
a radius of 2.85 fm. The calculations are carried out with the 
OLA. 
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FIG. 8: Prediction of the differential elastic-scattering cross 
section for p-^^O at _£;p= 300 MeV. The curve represents the 
full calculation with the harmonic-oscillator shell-model wave 
function. 



second minimum of the angular distribution. We have 
calculated the total reaction cross sections by varying 
the nuclear radii so as to be consistent with the observed 
interaction cross sections, and then analyzed the sensitiv- 
ity of the differential cross sections to the radii. We find 
a considerable change in the cross section, but whether 
it can really be observed or not strongly depends on the 
precision of the experimental data. 

We have predicted the differential cross section for p - 
elastic-scattering, which has recently been measured. 
Our prediction appears consistent with the preliminary 
data. This implies that it is possible to calculate to a 
good approximation the total reaction and elastic scat- 
tering cross sections at high energies. 
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APPENDIX A: 



HARMONIC-OSCILLATOR 
DENSITY 



We briefly explain the method of calculating the neu- 
tron orproton intrinsic density which is employed in our 
paper [3[. We assume a Slater determinant constructed 
from harmonic-oscillator wave functions. Denoting the 
Slater determinant by we first calculate the neutron 
or proton density which contains the effect of the center- 



of-mass motion 

A 

p(r)-(vI/|^<S(r,-r)P,|vI/), 



(Al) 



where is the nucleon single-particle coordinate and Pi 
is a projector to either neutron or proton. The single- 
particle wave function is given by the harmonic-oscillator 
function 



i^nljni = Rnlj{r)[Yi{r)xi/2] 



(A2) 



where Xi/2 denotes the spin function. If the orbit spec- 
ified by nlj is not completely filled, we take an average 
over the z component to, obtaining a spherical density. 
After some angular momentum algebra we have 



nlj 



(A3) 



where N„ij is the occupation number of the nlj orbit. 

To obtain the intrinsic density which has no contribu- 
tion from the center-of-mass motion, wc make use of the 
fact that the center-of-mass motion 4'cm(^) contained 
in 5* is factored out as 



* = «'0«'cm(X), 

where X is the center-of-mass coordinate and 



*cm(X) 



A 



3/4 



exp 



A 



(A4) 
(A5) 

(A6) 



Then the intrinsic density p(r) defined by 
A 

p{r) ^ {^^\Y,8{r,^ X - r)P,\^^) 
1=1 

is obtained from the following relation 

J dre^'^-^pir) = [(*cm|c"=-^|vl/„n)]"' 1 dr e*'= '^p(r). 

(A7) 

The radial functions which we need for the oxygen iso- 
topes are 



Roo{r) 
Roiir) 

Rw{r) 



1/2 

c 

1/2 



, 3 I °^PV O 2 



= ( 



16 \i/2 



r exp 



We give the examples of 21.23q ]-,giQ-^^^ neutron or 

proton density of ^^O is 
1 



p{r) 



2R^„ + 6i?^i + (M - 8)i?, 



in 



— ^-rexpi - — 



02 



1 + 2— 4 ^ ' 



15 



(A9) 
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were M is the number of protons Z=8 or the number of 
neutrons A^=13. The neutron density for ^'^O is 



2R 



■00 



■ cxp 



6R, 

7 

4 



02 



R 



10 



17 

15^ 



.(AlO) 



The corresponding intrinsic density is obtained from 
the relation (|A7p . The results are 



p{r) = 2 



3/2 



exp 



ur 



15 



1 + 3(1 - m) H a{l - u) 



-2M2(l + 5a(l-u)) 



for ^^O and 
p{r) = 2 



3/2 



exp 
17^3 



ur 

^2 



ri5 



- lOu 



17^2 



15 



(A12) 



3 3 

for 230. Here u = A/ {A - 1) and a = (M - 8)/15. 

As is evident, the above derivation assumes that the 
harmonic-oscillator constant v is the same for protons 
and neutrons. When we use different values for these, as 
in Table m the intrinsic densities of Eqs. (|A11|) and (|A12|) 
will no longer be exact. However, we have still used these 
equations to calculate the intrinsic densities in that case. 



APPENDIX B: ELASTIC SCATTERING 
AMPLITUDE 

In this appendix, following the method of Refs. we 
explain our calculation scheme of the proton-nucleus elas- 
tic scattering amplitude for a Slater determinant wave 
function constructed from harmonic-oscillator single- 
particle orbits. 

Let us start with the p -nucleus elastic-scattering am- 
plitude already given by Eq. ([1]), 

F{q) = ^ y" dbc'f^ {l - e*X"('')) , (Bl) 

where b is the impact parameter vector perpendicular to 
the beam (z) direction, K the initial momentum of the 
relative motion, and q the momentum transferred from 
the projectile to the target. 

As we have shown in Eq. ([2]), XN{b) is the nuclear part 
of the optical phase-shift function given by 



(*oin[ 



i-i^rp„(6-aO 



2 

1-^3. 



2 

(*o|0(b-s)|*o; 



l*o) 



(B2) 



where is the intrinsic (translation-invariant) A- 
nucleon wave function of the nuclear ground state. The 
coordinate Si is the projection onto the xy-plane of the 
nucleon coordinate relative to the center-of-mass of the 
nucleus, i.e., Si = s[ — X , where s'^ is the coordinate of a 
nucleon in the nucleus projected onto the aiy-plane, and 
X is the center-of-mass coordinate of the nucleus, ra. is 
1 for neutron and —1 for proton. For later convenience, 
here we introduce an operator 0{h — s). 

The main idea is that, for the harmonic-oscillator 
Slater determinant wave function, the elastic scattering 
amplitude F'(q) referring to the coordinate origin is fac- 
torized into a center-of-mass part and the intrinsic am- 
plitude F{q) 0, where F'{q) is defined by 



'(All) F\q) = ^ / rfbe"'-f'(*|5(3)(X-(l/^)yrj) 



27r 



x(l-0(b-s'))l*), 



(B3) 



For the harmonic oscillator, \^) is exactly factorized as 

|*> =|*c,n)|*0>. 

Then, the intrinsic amplitude which is free from the 
center-of-mass motion is easily calculated from F'{q) as 



F'iq) = (vI/em|e-'''-^|*cm>i^(g), 



(B4) 



where 



cxp[ 



-lyV/^A], (B5) 



and v is the size parameter of the harmonic-oscillator 
potential. The separability of the elastic scattering am- 
plitude, Eq. (|B3p . is valid when the size parameter of the 
harmonic well for protons, i^p, is equal to that for the neu- 
trons, i^n, which is the case for ^^O. On the other hand, 
when i/p j/„, which is the case for 20,21,22q^ only 
assume this separability and adopt the following relation 



i 



2 ^2 



N 2 



(B6) 



Now, we focus on the calculation of the elastic scat- 
tering amplitude F'{q) using the nucleon coordinates re- 
ferring to the coordinate origin. Since we represent the 
ground-state wave function as a Slater determinant. 



v[/^(A!)-i/2det[^™(r;)], 



(B7) 



the antisymmetrization is only required for the bra, such 
as 



(*! = (A!)-i/2det[5„,„^t^(r;)]. 



(B8) 



As one can see from Eq. (|B2[) . the operator 0{b — s') 
is factorizable into the operators that act on one particle 
subspaces, so that 

(*|0(6 - = det[0„™]„det[0„„]p, (B9) 



9 



where det[0„m]p is the Z y. Z determinant for Z protons 
(in this work, Z = 8) and det[0„TO]„ is the TV x TV de- 
terminant for TV neutrons. 0„m is the matrix element 
defined by 

0„,„ = <5„,„ - y V4(T-')r(b + s')^r^[r')dr' . (BIO) 
We have calculated the matrix clement, Onm, analyti- 



cally and obtained the determinants 0] ■ Concerning the 
neutron part, det[Onm]m we take the average of the de- 
terminants of different configurations. For example, for 
^''O, which has a 0+ ground state, we have 12 neutrons, 
four of which are in the 0^5/2 state. Therefore, we have 
three different determinants that give a 0+ ground state 
of 200. 
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